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INTRODUCTION 


oa  project  was  initiated  to  develop  a  methodology  to  predict  production 

Pt  ft  y  °  ^  giV6r  metal  Parts  line  through  computer  simulation.  At  the  pre¬ 

sent,  the  quantity  of  production  equipment  required  to  meet  a  specific  production 
rate  is  computed  on  the  basis  of  past  experience.  However,  it  is  apparent  that 

the^fi  lT  f°delS  f°r  S?Cl\  facilities  should  incorporate  today's  technology  in 
the  field  of  computer  simulation.  As  production  equipment  becomes  more  complex 

and  subsystems  are  added  to  the  facility  to  enhance  automation,  the  need  for 
computer  simulation  becomes  more  apparent  since  it  is  becoming  more  difficult  to 
pre  ct  production  capability  as  these  subsystems  are  added.  In  the  past,  a 
production  facility  consisted  mainly  of  semi-automatic  lathes  with  a  manual 

removal  P°Wered  conveyors,  manual  gaging,  and  manual  chip 

removal.  A  modern  facility  integrates  automated  machine  tools,  material 
handling,  chip  transport,  and  gaging. 

The  machine  tool  is  serviced  by  a  centralized  completely  automatic  handline 
system,  a  centralized  automatic  chip  transport  system,  and  automatic  gaging.  The 
subsystems,  material  handling,  chip  transport,  and  automatic  gaging,  can  directlv 
and  adversely  affect  production  throughput  and  the  monthly  production  capabilitv 
the  facility.  Past  methods  of  predicting  production  capability  counted  the 
number  of  machine  tools  per  operation  and  related  production  rates  to  forecast 

f acilitv°n  f  Vhole*  This  can  ^nger  be  done  in  a  modernized 

facility.  With  the  introduction  of  automated  systems  into  a  facility,  personnel 

are  no  longer  relied  upon  to  maintain  throughput  by  manually  moving  shell 
shoveling  chips  or  manually  gaging  parts.  Emphasis  is  now  placed  on  these  auto¬ 
mated  systems  and  must  be  considered  in  predicting  throughput  and  monthly  produc- 
tion  capabilities  of  modern  facilities. 

This  report  covers  the  simulation  of  various  metal  parts  lines  usine  exist¬ 
ing  simulation  models  as  well  as  the  development  of  alternate  models  with  which 
various  line  characteristics  and  configurations  are  developed  and  tested.  The 

“  °f  aCqUhirl"g,  the  data  for  simulation  that  was  not  readily 

sllvP^  ff  /  7  18  addressed  as  wel1  a*  developing  a  methodology  to 

solve  buffer  sizes  for  automated  material  handling  systems  necessary  if  any  accu¬ 
racy  is  to  be  obtained  from  the  prediction  model.  7 

n„JhrrettCal  resalts  obtained  from  these  various  simulation  models  are  com¬ 
pared  to  actual  production  rates  and  design  rates  to  which  several  conclusions 
and  recommendations  will  be  discussed. 


INVESTIGATION  OF  VARIOUS  SIMULATION  METHODS 


Simulation  has  been  defined  as  the  process  of  designing  a  model  of  a  real 
system  and  conducting  experiments  with  this  model  for  the  purpose  either  of 
understanding  the  behavior  of  the  system  or  of  evaluating  various  stratiiM 

of \h“Sy8t.M  rred  by  a  crlterl0n  0r  set  o£  criteria)  for  the  operation 
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were  I^ead^ivil^lT  thlS  /'““.“l  ““  to  deter"i‘“  what  simulation  models 
!rl!a  f  ^  available.  It  was  found  that  a  profusion  exists  in  simulation  pro- 

g  ms  for  production  lines.  The  problem  was  in  finding  and  selecting  the  proper 

program  for  our  particular  application.  Two  programs  appeared  to  be  most  suited 

System  r(Gpls)S>  and6  th3S  th  SimfUlati°n  lanSuaSe  titled  General  Purpose  Simulation 
ab!a2LiGPS»  u  d  th  0ther  General  Modeling  (GENMOD)  developed  and  in  use  at 
ARRADCOM.  Both  GENMOD  and  GPSS  are  simulation  models  that  use  software  packages 

f  ™'ha  USeVS  W—  —  provide  various  output  iTZ- 

Mred  to  oLi„»,  generad  purpose  simulation  models  and  languages  which,  when  def¬ 
lation  langu^fnc/”^^  p™8rar"ln8  languages,  differ  In  that  the  slmu- 
occur*1  1S  TT  ,T  controlling  the  sequence  in  which  events 

provide  mh  ? 1Vantage  of  the  two  models  is  that  the  simulation  languages 

capacity  violating?8”0  Pr°8ramS  “hlCh  ^  £°r  l0glCal  e"°r'  ayta*'  aad 

r^SimUlatiu°n  models  may  be  either  continuous  change  or  discrete-event  tvoes 
a  continuous^ flow  7 i %  appropriate  in  simulating  a  system  which  consists  of 
in  individual  terms  iTfOT^tioa  or  niaterial  counted  in  the  aggregate  rather  than 
tualized  Is  ;  discrete-event  models  changes  in  the  system  are  concep¬ 

tualized  as  distinct  occurrences;  parts  flow  through  the  system  requiring  the 

?LlOrma0nuCeeuitathe°PerathiOn  °n  Part  bef°re  U  Can  —  °n  t0  the  opera- 

tion.  Queuing  theory  becomes  important  in  that  parts  may  wait  in  line  before 

thonrcsuuiCei„dicIntf0mha';10i;  regardlng  the  behavlor  of  the  model  Is  tabulated  and 
GPSS  and  rFmnn  a  •  hapPens  under  a  Particular  set  of  assumptions.  Both 
appropriate  ^han316  dascrete-d''e“  <*pe  simulation  models  and  are  considered  more 
lines  Continuous  change  models  when  simulating  metal  parts  ammunition 


able  at  °”. ,  1"’'estl8a“on  of  the  various  simulation  programs  found  avail¬ 
able  at  the  time,  the  GENMOD  computer  program  was  selected  as  the  most  appro- 

priate  program  for  our  intended  use.  It  is  a  generalized  model  which  can  inves¬ 
tigate  various  automated  production  lines  and  the  efforts  of  proposed  design 
changes  without  constructing  an  entirely  new  model,  thereby  saving  costly  time  in 
making  models  unique  to  our  production  line  alone  (ref  2). 


DATA  ACQUISITION 


rafo  aay,  simulation  Program  to  be  accurate,  the  data  input  must  also  be  accu- 

nroiect  Thl  PreS*nted  °ne  of  the  maJor  Problems  at  the  very  beginning  of  this 

available  M  f  data  re(*uired  for  most  of  the  simulation  programs  was  not 
available.  Most  simulation  programs  under  consideration  for  our  application 

theUtred  maan'tlme-between-fail^re  (MTBF)  and  mean-time-to-repair  (MTTR)  data  for 

^cuLrlv^  eqUlpmant  tD  be  simulated*  A  review  of  data  in  the  industry  par- 
,lyf  1  Projectile  metal  parts  plants  revealed  such  data,  at  least  in  the 

fhb.re  ?rm>  W3k  n0t  available*  Most  plant  operators  knew  the  availability  of 
their  equipment  but  not  the  MTBF  or  MTTR.  Availability  is  defined  here  as  the 

:  Un±t  °f  productioa  equipment  will  operate  when  re- 
able  o^MTBF  pl„sCamTR  dafl”ad  “  ^  dlvldad  by  tha  tdtal  C1“a  a''ail- 
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Availability  = 


MTBF 

MTBF  +  MTTR 


Although  fairly  reliable  data  on  availability  of  equipment  based  on  past 
experience  could  be  readily  obtained,  the  MTBF  and  MTTR  needed  to  run  any  simu¬ 
lation  program  could  not  be  accumulated.  It  was  apparent  that  MTBF  and  MTTR  data 
had  to  be  collected. 


All  published  reports  were  screened  to  determine  the  available  data.  Masor 
Chamberlain  Incorporated,  the  operator  of  a  new  facility  now  under  construction 
had  contracted  with  GARD  Incorporated  for  a  report  (Sept  78)  to  the  Project 
Manager  Base  Modernization  office  to  analyze  the  Mississippi  Army  Ammunition 
Plant  (MSAAP)  complex  using  the  GENMOD  program.  GARD  encountered  the  same  pro¬ 
blem  of  attaining  MTBF  and  MTTR  data.  GARD  reviewed  downtime  records  from  the 
various  Chamberlain  plants  at  Scranton,  Pennsylvania,  New  Bedford,  Massachusetts 
and  Waterloo,  Iowa  and  developed  MTBF  and  MTTR  values  (table  1)  for  the  MSAAP 
faciiity  for  the  M483  projectile  metal  parts.  GARD  then  ran  the  GENMOD  program 
and  the  results  were  encouraging  in  the  respect  that  the  design  rate  for  the  line 
was  predicted  using  reasonable  buffers  on  the  automatic  material  handling  system 
of  320  spaces  in  all  but  the  forge,  heat  treat,  and  phosphate  and  lube  banks. 


A  review  of  the  data  together  with  the  associated  availability  leads  to  some 
questions  as  to  the  reliability  of  the  data.  All  of  the  availability  numbers 
appear  to  be  too  high  based  on  past  experience.  For  instance,  all  availabilities 
except  for  forge  equipment  are  around  95%,  whereas  the  expected  availabilities 
are  in  the  order  of  80%  or  possibly  lower  in  some  cases. 


Because  of  this  doubt,  it  became  apparent  that  more  information  with  greater 
accuracy  was  required  for  an  accurate  simulation  to  be  achieved.  Fortunately 
the  National  Presto  Industries  (NPI)  in  Eau  Claire,  Wisconsin  that  produced  the 
105-mm  Ml  HE  projectile  in  a  large  volume  was  known  to  keep  good  maintenance 
records.  The  company  was  visited  with  the  intent  of  obtaining  equipment  oper¬ 
ating  records  for  statistical  analysis  yielding  MTBF  and  MTTR  values.  Under 
government  auspices,  NPI  conducted  a  rebuild  and  modernization  program  on  their 
number  one  forge  line  which  culminated  in  an  equipment  proveout  in  the  spring  of 
1973.  During  this  proveout,  extensive  downtime  records  were  taken  extending  over 
a  period  of  2  months  during  which  the  line  ran  24  hours  a  day,  5  days  a  week.  An 
example  of  the  record  sheet  is  shown  in  table  2.  Information  of  this  caliber 

proves  invaluable  as  far  as  input  data  necessary  to  running  a  computer  Simula- 
tion. 


.  fo.rge  line  consisted  of  six  operations  in  series  with  each  beginning 
with  induction  heating  the  billet  to  approximately  1060°C  (2100°F).  The  billet 
then  goes  through  a  high  pressure  water  descaler  to  remove  iron  oxide  buildup  on 
the  surface  caused  by  heating.  The  block,  cabbage,  and  pierce  forging  operations 
on  one  transfer  press  are  than  performed  on  the  billet  producing  a  bottle  type 
configuration.  The  forged  shell  is  then  control  cooled  to  room  temperature  to 
produce  uniform  metallurgical  properties  throughout  the  shell. 

A  failure  event  in  this  line  was  considered  to  be  any  repair  needed  in  any 
one  of  the  six  operations.  This  failure  results  in  a  loss  of  production.  Var¬ 
ious  uptimes  (MTBF)  and  downtimes  (MTTR)  were  calculated  from  the  first  two  col¬ 
umns  of  table  2.  This  forge  ran  24  hours  without  any  stops  for  lunch  or  breaks; 
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therefore,  any  stop  in  the  line  was  due  to  failure  in  machine  operation.  To 
illustrate  a  calculation,  the  table  shows  that  at  4:00  the  forge  went  down  and 
started  up  at  4:01  producing  a  downtime  event  of  1  minute,  then  ran  from  4:01  to 
4:35  producing  an  uptime  event  of  34  minutes.  The  same  analysis  was  used 
throughout  the  study  covering  two  months  in  which  approximately  1600  individual 
uptime  and  downtime  events  were  recorded. 

The  time-to-failure  and  time-to-repair  events  were  totaled  and  categorized 
into  time  intervals  for  graphical  purposes.  By  individually  graphing  the  fre¬ 
quency  of  repairs  versus  the  time-to-repair,  histograms  were  produced  which  can 
visually  be  interpreted  to  follow  negative  exponented  functions. 

The  J-fi.r1st  histogram  (fig.  1)  is  a  graph  of  frequency  of  failures  versus  time 
between  failures  which  shows  that  approximately  180  times  the  forge  line  was 
running  between  0  and  15  minutes  before  a  repair  was  needed.  The  next  interval 
shows  that  approximately  150  times  the  forge  was  running  between  16  and  30  min¬ 
utes  before  a  repair  was  needed,  etc.  The  overall  pattern  appears  to  be  a  nega¬ 
tive  exponential  function.  _  The  MTBF  was  calculated  to  be  59.08  minutes  and  is 
defined  as  average  MTBF  =  X  =  -,  where  X  =  individual  MTBFs  and  N  =  number  of 
occurrences* 

The  second  histogram  (fig.  2)  is  a  graph  of  frequency  of  repairs  versus  time 
to  repair  which  shows  that  around  350  repairs  were  needed  taking  between  0  and  5 
minutes.  The  second  interval  decreases  sharply  in  number  of  repairs,  150  repairs 
taking  between  6  and  10  minutes.  Again  the  pattern  appears  to  be  a  negative 
exponential  function.  The  MTTR  was  calculated  to  be  13.64  minutes. 

These  figures  do  not  include  catastrophic  failures  which  alter  the  two  fig¬ 
ures  to  53.02  min  and  27.27  min  for  MTBF  and  MTTR,  respectively.  The  forge  line 
runs  at  an  availability  of  82%  until  a  catastrophic  failure  occurs  which  then 
drops  the  availability  to  approximately  65%  overall.  A  catastrophic  failure  is 
defined  as  a  failure  which  causes  a  loss  of  production  whose  magnitude  extends 
significantly  beyond  the  range  of  normal  expected  values.  Experience  has  taught 
production  and  manufacturing  engineers  in  these  plants  to  bank  one  week  of  pro¬ 
duction  forgings  in  order  to  supply  machines  down  the  line  with  parts  to  compen¬ 
sate  for  the  catastrophic  failures  while  repairs  are  being  made  on  the  faulty 
piece  of  equipment. 

/ 

The  theory  that  both  histograms  are  exponential  in  nature  was  tested  statis¬ 
tically  by  using  the  Kolmogorox-Smirnov  Goodness  of  Fit  Test  (table  3).  This 
test  is  based  on  a  comparison  of  the  cumulative  frequency  distribution  specified 
by  the  theory  to  the  observed  cumulative  frequency  distribution.  Therefore,  this 
test  compares,  interval  by  interval,  the  areas  of  a  theoretical  histogram  gener¬ 
ated  by  the  specified  function  to  the  observed  histogram.  The  test  states  that 
when  subtracting  the  calculated  from  the  observed  a  good  correlation  exists  if 
the  difference  is  close  to  zero. 

Dx  =  max  [F(x)  calc  -  F(x)  obs]  (ref  3). 

All  values  approach  zero  indicating  that  distribution  is  close  to  the  assumed 
exponential  (table  3). 
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Test  results  show  that  the  failure  cycle  and  the  characteristics  of  the 
forge  line  can  be  identified  by  the  negative  exponential  function  of  the  form 

„  *  “At  ^  1 

Y  =  Ae  where  A'x -  and  e  =  epsilon 


mean 


Time-to-f ailure  equation 

-t/ 59.08 


Y  = 


1 


59.08 

Time-to-repair  equation 


-t/ 13.64 


These  equations  are  important  for  computer  simulation,  since  events  occurring  in 
time  have  now  been  shown  to  have  a  pattern  of  occurrence  which  can  be  imitated, 
thereby  increasing  knowledge  of  line  operation  and  design.  However,  in  order  to 
simulate  a  line,  added  information  is  needed  on  equipment  other  than  just  the 
forge  line ;  i.e.,  machine  tool  turning  equipment. 

The  opportunity  to  collect  data  on  machine  tools  came  during  the  proveout  of 
the^  155-mm  M483  metal-parts  line  at  Norris  Industries  and  was  collected  for  a 
period  of  20  days,  8  hours  a  day. 

Data  sheets  used  during  the  proveout  (table  4)  are  similar  to  those  used 
during  the  forge  line  proveout,  except  for  the  column  marked  event  code.  This 
column  is  used  by  the  data  collector  to  identify  what  type  of  event  caused  loss 
of  production  to  the  machine  and  is  broken  down  into  eight  categories  as  can  be 
seen  on  the  bottom  portion  of  the  sheet. 


A  computer  program,  developed  by  ARRADCOM's  Quality  Assurance  Division, 
analyzed  the  proveout  data  and  provided  information  detailing  machine  breakdowns * 
tool  changes,  and  observed  production  rates.  An  example  of  a  computer  printout 
of  the  daily  downtime  record  sheets  is  shown  as  tables  5  and  6  with  table  7  pro¬ 
viding  a  summary  for  a  particular  machine  for  the  overall  test.  A  system  summary 
is  also  provided  (table  8)  analyzing  by  machine  and  by  operation  various  machine 
operating  characteristics. 


.  lWith  the  acquisition  of  this  data  (MTBF  and  MTTR  values)  is  was  possible  to 
simulate  metal  parts  ammunition  lines  by  computer. 


SIMULATION  OF  VARIOUS  AMMUNITION  LINES 


Norris  Industries 


The  first  line  to  be  simulated  was  the  Norris  Rough  and  Finish  Turn  M483 
line  whose  input  data  is  listed  in  table  8. 
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operations Tade  Unlimited1  (7e  “STt  f S  ^  C°ntr01  P°lntS  <bu£fa->  batwee" 

rt:  :^h^5;!>?r  *s=  “= 

specif  icaUjT  ^STm  ~  ^  “  <“b“  ^SSJLfSS  Te 

sBXsS  ESHF“  -s  j^.: 

SlsSrSoSffet  ~  »°^rinl8thtf  c"asde  T  BII  e^ls  ^  322?“, h°e 

releases  11.7  psrts/BTI  (46!s  parts/hr)'  ihe^e"8  ?”  BTI  U0°-4  parts/hr)  a"“ 
line  is  at  fault  'i  V7  parts/hr).  The  decision  as  to  which  side  of  the 

case  is  is  nnn  ,  fS  t0  thG  intended  d^ign  rate  of  the  line,  which  in  this 
nartc  15»000  acceptable  units  per  500-hour  month.  An  average  rate  of  30  net 

buffer  ?in  ea^h  ITse  the  Pr°blem  *S being  an  overfed 

Turn  ID-OD’)  anH  un  rv  *  atlon  of  one  machine  m  both  operations  210  (Rough 

Auction  \.t?  since  ‘  th"e  S  "““K  “m  not  *f£e<:t  ‘he  ‘"tended  designed 

tSaterSatS^ 

its  Srt*:  SSs/SSs'^ - 1 « 


Projected  average  outputs 

per  basic  interval 
per  27  interval  shift 
per  3/8/5  week 
per  63  shift  month 


9.86  (0.66  per  min) 
270.0 
4043.0 
16979.0 


^rzl&Ss&pii-ZsZ'  ««=  «=s 

tomSSTrSs  :fLSi;ss  if  rrs  ln  che  une 

ten7VidUedtUups1t°rI‘eamhebef£*re  con8te«ly billing  upl  shoos'^hatT’oVmaUy  lo»“con- 

~  :rrs7;ss7e- 

the  welT  heat  St,’  SZZoTJtLT  “  ha"dle  ^  b“ff»= 
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Another  way  of  analyzing  the  projected  rates  is  to  compare  it  with  actual 
production  rates  for  accuracy  within  the  simulation  run.  This  can  be  done  bv 
comparing  production  rate  of  operation  305,  the  last  operation  in  the  simulation 
to  the  proveout.  By  using  the  proveouts  daily  summary  sheets  (tables  5  and  6) 
and  relating  total  parts  processed  to  scheduled  uptime,  an  overall  rate  of  37  95 
parts  per  hour  results.  GENMOD  identified  machines  36  and  37  in  the  simulation 
to  be  those  of  operation  305,  and  by  calculating  the  total  parts  processed  by 
each  machine  to  the  scheduled  production  run,  a  rate  of  37.88  parts  per  lour 
results.  Both  rates  correlate  well  which  lends  to  a  high  degree  of  confidence  in 
the  output  information  generated  by  GENMOD  on  the  Norris  155-mm  M483  line. 


Mississippi  Army  Ammunition  Plant  (MSAAP)— Mason  Chamberlain  Inc. 


The  next  facility  to  be  simulated  was  the  MSAAP  155-mm  M483  line  (fie  31 
currently  under  construction  and  operated  by  Mason  Chamberlain,  Inc.  Production 
rates  were  obtained  from  equipment  proposals  while  MTTR  and  MTBF  values  w=>re 
extracted  from  similar  equipment  at  both  Norris  and  National  Presto  (table  12).~ 

to  Si“ulat,i°n  °*  th®  line  allowed  all  buffers  between  each  operation 
to  be  infinite.  Downtimes  for  large  equipment  such  as  stress  relieve  and  cool 

125S^aCC  anC  a  ’  anC.clean  ”ere  Set  at  1  hour  with  uptimes  ranging  from  160  to 
,  u  *  y  01ng  this  the  line  would  be  run  at  peak  production  output  without 

otherVff  la!rSe,  equipment  beinS  down‘  This  Places  a  greater  strain  on 
nl  7n  7  P  S  n°  limit  uP°n  each>  imbalances  in  the  line  will  show 
up  in  the  buffer  content  summary  as  large  numbers.  A  list  of  the  buffers  and 
operations  which  surround  each  are  shown  in  table  13. 


,  ®uffers  4,  8,  and  15  are  located  around  large  equipment  such  as  phosphate 
and  clean.  Here  each  buffer  built  up  its  contents  even  though  the  large  machin- 

Ceed  e?PerKenw  Sh  a™ilabilities  with  low  downtimes.  This  illustrates  the 

need  for  banking  at  such  areas  in  metal  parts  lines  which  is  exactly  what  facil- 
ity  design  engineers  do  based  upon  experience. 

Buffer  2  which  is  between  operations  shot  blast  and  concentric  turn  has 

is  ch  1  I"'  UPs  <U  throuShout  the  simulation.  Part  of  the  reason  for  tUs 

fLl/hrl  gTnPr°,nC  “  rate  °f  Sh0t  blaSt  <396/hr>  versus  concentric  turns 

Idol/ hr ).  The  other  reason  is  concentric  turn’s  MTBF  is  half  that  of  the  slot 

blast  value  and  therefore  experiences  a  buildup  of  parts  in  buffer  2. 

l  AoaBetWten  °^tti0nS  b°re  n°Se  and  rouSh  turn  body  is  buffer  7  which  averaged 

i^bal  an  '  A  a  ^  “  Valu6S  f°r  each  °Perati™  were  made  the  same  so  ? he 
imbalance  must  be  due  to  production  rate  differences  of  360  per  hour  versus  324 
per*  tiour . 


u  Bufife^  21  ^ls°  haS  a  high  averase  of  211  parts.  Both  operations  (Fiberglass 

MTBF  values11  Tt  ^ibe^glass)  hava  the  same  production  rates  and  in  examining  the 
MTBF  values  it  is  shown  that  the  wrap  machines  only  stay  running  half  that  of  tne 

turn  equipment.  Why  does  buffer  21  fill  up?  The  answer  can  be  found  by  examin- 

(Turn  Fib  \  T  failures  for  each  operation  are  listed.  Operation  1155 

(Turn  Fiberglass)  experienced  a  total  of  579  failures  for  an  availability  of 
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83.93%.  Operation  1160  experienced  464  failures  for  an  availability  of  88  827 

ITT  ;  U,65  had  m°re  failureS’  buffer  21  *ad  build  up 1  i ts  content 
The  quoted  production  rate  for  the  simulation  was  138,908.  contents. 

r;;„-  r-Sr;*! 

Jieldiiri26rf84eCptaSrtsUP°n  500  hPr°jeCted  °Utputs  with" Ihe^l-Cu^dov^t ime^ run 
only  103  884  „*  L*  PfL 51?0_hour  month»  while  the  8-hour  downtime  run  projects 

P^o™^8.  per  500‘hour  “nth-  "»  desi*"  ““ 

buffed  /rrav^Lu%air^%rreeS°;e-the  Msr  .iine  usi^  “«•»« 

table  15  A  vr^nh  nf  IT  values  with  the  resulting  production  rates  is  shown  in 

buffer  size  between  operations  (  Th  ^  Sh°WS  production  rates  as  a  function  of 
equipment  at  a  m/nTm  0peratl°ns*  The  importance  in  keeping  downtimes  of  large 
quipment  at  a  minimum  is  shown  in  runs  1  and  3  (table  15)  which  keep  the  total 

to  8fcWsihrproieCc0tnSdtan,:  H(130)  ^  °peration  and  varies  downtimes  from  1  hour 
nours  the  projected  production  rate  drop  from  127  K  to  105  K. 

figurre4/efmeteiargereltahe°bufXrS  producrtlon  ,tat?  and  buffer  sire  Is  shown  In 

8  hours  ™  £ 
buffers  to  be  as  low  as°  lT.nT.UU  ‘  h°Ur  “  ^  all<”' 

livens 


Louisiana  Army  Ammunition  Plant — Thiokol  Corp. 


£g-TlF  fin  ““ 

sizes.  X^LuicTbe^th  “T?  dlfferent  MTBF  and  MTTR  values  and  buffer 

unlimited  buffers  were  usedic^S^  Xh"/  computar  simulations  using  GENMOD, 
the  line  In  rh P  K  Jl  f  1  identify  aay  imbalances  or  critical  areas  in 

q  aii  ,  ,  average  buffer  content  summary  (table  17),  buffers  1,  3  7  and 

frr, 

i--ot-  —  -  -  -  -  c~ 


8 


Runs  were  made  varying  the  MTBF  and  MTTR  values  (table  18)  on  operation  1110 
(slow  cool)  and  on  operation  1270  (heat  treat).  As  predicted,  the  runs  with  long 
repair  rates  yielded  lower  monthly  production  rates  when  compared  with  runs  which 
used  short  repair  rates.  This  is  true  even  when  the  availabilities  are  made 
higher  for  long  repair  rates  when  compared  with  short  repair  rates  and  lower 
availabilities  on  operations  such  as  heat  treat  and  slow  cool.  An  explanation  is 
that  key  operations  (1  machine)  such  as  slow  cool  or  heat  treat  have  a  great 
affect  on  production  when  experiencing  long  repairs  causing  depleted  downstream 
inventories  and  requiring  extra  time  in  replenishing  buffers  to  stable  capaci- 
ties. 


Buffer  limits  also  place  restraints  on  predicted  production  rates  to  an  even 
greater  degree  than  that  of  MTBF  and  MTTR  values  (table  19).  By  placing  limits 
of  500  (at  forge  and  heat  treat  areas)  and  88  (actual  buffer  sizes)  on  the  line 
the  production  rates  falls  from  59,717  to  51,384.  However,  run  BAESF80  used  the 
large  downtime  figures  of  158  minutes.  By  changing  this  number  to  56  minutes 
(run  BAESF32)  the  production  rate  can  be  raised  to  58,581  per  month  with  the  same 
buffer  limit  of  500  and  88. 

These  computer  simulations,  show  the  various  f kctuatioas  in  production 
levels  that  occur  based  upon  ^nging  uptimes  and  downtimes  as  well  as  buffer 
sizes.  The  computer  simulations  of  the  LAAP  line  show  that  downtimes  of  the 
magnitude  of  189  minutes  (on  key  operations)  will  prevent  achieving  the  designed 
production  rate  of  63  K  per  month,  but  if  reduced  to  approximately  60  minutes, 
production  approaches  the  design  rate.  This  is  true  even  if  the  uptimes  are  as 
large  at  790  minutes  and  a  low  as  158  minutes.  By  decreasing  the  downtime  of 
such  operations,  buffer  sizes  of  500  and  88  will  work  effectively. 


MATERIAL  HANDLING  BUFFER  SIZE 


Along  with  MTBF  and  MTTR,  buffer  size  must  also  be  provided  as  input  to  any 
simulation  program.  The  program  selected  for  our  studies,  GENMOD,  is  no  differ¬ 
ent  and  a  methodology  to  derive  the  proper  buffer  size  became  highly  desirable. 
Initial  efforts  in  establishing  realistic  buffer  sizes  consisted  of  contacting 
material  handling  system  contractors.  These  efforts  revealed  that  no  system 
exists  and  that  buffer  sizes  are  established  rather  arbitrarily  by  the  buyer  of 
the  material  handling  system  (MHS)  and  in  some  cases  by  the  seller  of  the  MHS 
Since  there  has  been  very  little  experience  in  the  projectile  metal  parts  busi¬ 
ness  with  automated  MHS,  review  of  past  histories  in  existing  facilities  was  not 
very  helpful.  Therefore,  initially  the  simulation  programs  were  run  based  on 
either  infinite  buffer  sizes  or  some  arbitrary  buffer  size  based  on  cost  conside- 
ration  and  very  limited  experience. 


If  any  simulation  program  was  to  be  reliable 
developed  to  derive  buffer  sizes  based  on  some 
methodology  developed  and  used  for  this  program  is 


,  a  methodolgy  would  have  to  be 
mathematical  principles.  The 
presented  here. 


M  **  j£°Uld,  hl  understood  what  the  function  of  a  buffer  is  and  how  it  relates 
to  the  MHS  and  the  plant  operations.  It  should  be  noted  that  when  we  have  sev¬ 
eral  machines  m  each  operation  these  machines  are  turning  off  and  on  randomly. 
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based  on  their  availability.  The  real  time  production  rate  at  each  operation  in 
a  sequence  of  operations  varies  accordingly.  In  such  a  situation  the  MHS  must 
have  buffer  to  meet  one  of  two  conditions.  Consider  the  case  where  you  have  a 
MHS  between  operation  A  and  operation  B. 

Operation  Operation 


Buffer 


A 


If  Operation  A  for  a  short  period  of  time  is  producing  more  parts  than  B,  the 
buffer  (MHS)  must  absorb  parts  from  operation  A  until  operation  B  can  accept 
these  parts.  On  the  other  hand  if  operation  A  for  a  short  period  of  time  is 
producing  less  parts  than  B,  the  buffer  (MHS)  must  provide  parts  to  operation  B 
until  operation  A  comes  back  up  to  speed.  These  variations  (speed  or  rate)  of 
output  of  A  and  B  are  caused  by  the  random  turning  on  and  off  of  several  machines 
in  each  operation  in  accordance  with  their  respective  availability.  It  was 
stated  earlier  that  availability  which  is  more  readily  known  for  production 
equipment  consist  of 


Availability  = 


MTBF 

MTBF  +  MTTR 


Unlike  the  situation  for  GENMOD  where  we  must  know  the  actual  values  for 
MTBF  and  MTTR  for  each  type  of  equipment,  the  availability  (A)  is  sufficient 
information  for  calculating  buffer  size  for  the  methodology  presented  here.  If 
it  is  stated  that  a  machine  tool  has  an  availability  of  80%,  it  is  expected  to 
operate  80%  of  the  time.  This  is  a  life  or  death,  yes  or  no  situation.  A 
machine  is  either  running  or  it  is  not  and  no  other  condition  exists.  Further, 
these  two  events  are  exclusive  of  each  other  and  together  these  two  events,  on  or 
off,  make  up  the  complete  time;  no  other  event  exists.  Such  a  situation  lends 
itself  to  the  binomial  distribution  (fig.  6)  which  was  used  to  develop  buffer 
sizes  for  a  specific  application  or  production  line.  The  following  explanation 
is  an  example  of  the  precedure  followed. 

Consider  the  case  where  you  have  a  line  consisting  of  a  series  of  operations 
where  you  want  a  production  capability  of  65,000  units  per  500  working-hour 
month.  Therefore,  the  desired  throughput  at  any  time  would  be  65,000  divided  by 
500  or  130  units  per  hour. 


65,000  units/month  , ^ 

500  hour/month  =  130  units/hour 

The  130  units  per  hour  is  the  mean  throughput  value. 

Now  consider  the  situation  where  there  are  two  operations,  A  and  B.  In 
operation  A  assume  five  machine  tools  are  required  to  meet  the  production  rate 
and  six  machine  tools  are  required  in  operation  B.  In  addition,  assume  that  the 
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availability  of  these  machine  tools  is  80%  or  0.8.  Therefore,  in  operation  A  (5 

machines)  26  parts  per  hour  net  or  32.5  parts  per  hour  gross  must  be  produced  for 
each  machine. 


or 


130  parts/hour 

“5  machines  =  26  P*rts/hour /machine  (net) 


26 

0.8 


32.5  parts /hour/machine 


The  actual  cycle  time  for  the  machine  in  operation  A  is  110.77  seconds/part . 


3600  seconds/hour 

"32. 5  parts /hour  =  110,77  seconds/part 


In  operation  B  (6  machines)  21.667  parts  per  hour  net  or  27.0833  parts  per 
hour  gross  must  be  produced  for  each  machine. 


130  parts/hour 

~6  machines  =  21 ,667  Parts/hour /machine  (net) 


21.667  _ _ 

Q#g  ”  27.0833  parts/hour/machine  (gptfss) 

The  actual  cycle  time  for  the  machine  in  operation  B  ^  132.92  seconds/part. 
3600  seconds/hour 

27.0833  parts/hour  =  132,92  seconds/part 


these  machines  randomly  turn  on  and  off  in  each  operation  but  at  an 
availability  of  80%,  the  percent  of  time  a  specific  number  of  machines  is  oper- 
ating  can  be  determined  by  using  the  binomial  distribution.  For  operation  A 
(table  20)  exactly  five  machines  will  be  operating  32.768%  of  the  time,  four 
machines  will  be  operating  40.96%  of  the  time,  or  at  least  four  machines  will  be 
operating  73.728%  of  the  time.  The  expected  time  for  the  remainder  of  the 
machines  can  be  calculated  by  using  the  binomial  distribution.  The  time  for 
operation  B  with  six  machines  can  be  calculated  by  the  same  method  (table  21). 


By  taking  a  unit  of  time  (hr)  the  various  rates  of  production  over  the  per¬ 
centage  of  time  the  machines  are  operating  (up)  can  be  calculated  for  that 
°ur*  F°r  axamPle>  in  operation  A  five  machines  are  up  and  running  32.768%  of 
the  time;  therefore,  parts  are  being  produced  at  the  rate  of  162.5  parts  per 
hour.  This  occurs  over  a^fime  period  of  19.6608  minutes.  / 
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(5  machines)  (32.5  £----?/hour  = 

machine 


162.5  parts  per  hour) 


(0.32768)  (60  minutes)  —  19.6608  minutes 


Since  162.5  parts  per  hour  is  2.70833  parts  per  minute,  operating  for  19.6608 
minutes  produces  53.248  parts. 


no  c* c* aq  _  (162.5  parts/hour) 

(.19.0008  minutes  — 77-7 - - - ■ - L  -  co  040 

(60  minutes/hour)  parts 

The  above  calculations  are  shown  in  tables  22  and  23,  columns  I  through  V.  These 
co  umns  can  be  completed  by  using  similar  calculations  for  each  row  developed  for 
opera  ons  an  B  starting  from  the  top  row  and  working  down.  Columns  VI 

through  X  for  each  operation  can  be  calculated  by  starting  from  the  bottom  row 
ana  working  up. 

1...  Tt  1S  unllkely  that  machines  will  follow  either  pattern,  but  these  two  con- 
!  h  n  fh?MSen!  the  extreme  of  all  other  conditions  that  will  occur  randomly 
within  the  binomial  distribution.  A  definition  of  a  maximum  condition  is  when 

!  "  Hg  °T  With  311  machines  in  a  given  operation  up  and  running  and  a  minimum 
condition  is  when  starting  out  with  none  of  the  machines  operating.  Therefore 
columns  I  through  V  represent  the  maximum  condition  where  the  majority  of  parts 
are  pro  uced  at  the  front  end  of  the  time  interval  and  columns  VI  through  X 

f  rh  ml.n/mU“  condition  where  ^e  majority  of  parts  are  produced  at  the 
back  end  of  the  time  interval. 

Operation  A  producing  at  the  maximum  and  operation  B  producing  at  the  mini¬ 
mum  are  plotted  in  figure  7.  The  figure  shows  that  at  the  time  interval  Tto 
rh  S  °paration  A  Produces  more  parts  than  operation  B  can  accept;  therefore, 
c  between  the  two  operations  must  be  able  to  accept  or  buffer  these 
parts.  The  maximum  vertical  distance  between  the  plot  of  the  two  operations 
represents  that  buffer  and  measures  approximately  20  parts.  Since  the  maximum 

and  minimum  conditions  represent  the^Wfemes  of  all  the  situations  that  might 
occur,  a  buffer  size  of  less  than  20  parts  absorption  is  actually  required  in 

Tn^MT^R8*  question  arises  as  to  what  is  the  expected  frequency  that  MAX  A 

and  MIN  B  would  occur.  MAX  A  must  have  a  special  case  of  5,  4,  3,  2,  1  machines 

up  equ  va  ent  or  similar  to  a  5— card  draw  without  replacement  equal  to  5!.  MIN  B 

must  ave  0,  1,  2,  3,  5,  6  equal  to  a  6-card  draw  without  replacement  or  6!.  The 
requency  of  a  MAX  A  and  MIN  B  occurring  simultaneously  would  be  the  product  of 

ft^Ann0  r  (  °  (6°  =  86’40°*  Therefore ,  this  event  could  occur  once  every 

,  t  me  intervals.  Since  the  time  interval  is  1  hour,  a  MAX  A-MIN  B  could 
occur  once  every  86,400  hours.  Some  cases  may  not  have  five  and  six  machines 

machinl^ofvor’N  _0perationf  such  as  forge  or  heat  treatment  may  have  only  two 
[(2.) (2.)  4  hours] .  In  this  case  the  maximum  or  minimum  condition  are 

met  every  4  hours.  If  a  40-hour  capability  without  exceeding  buffer  capability 
is  needed,  take  the  buffer  quantity  calculated  for  one  hour  time  interval  and 

multiply  that  value  by  40  hours  divided  by  4  equals  10.  In  many  instances  this 
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low  number  of  machines  coupled  with  the  complexity  of  the  machines  resulting  in 
availabilities  lower  than  0.8  has  led  to  manufacturers  banking  as  much  as  one 
week's  production  of  process  pieces  to  maintian  throughput  to  meet  monthly  pro¬ 
duction  deliveries. 


Consider  again  operation  A  and  B  where  the  opposite  of  the  above  example 
occurs.  Consider  the  situation  where  operation  A  produces  at  the  minimum  and 
operation  B  produces  at  the  maximum. 


A  similar  situation  exists  for  MIN  A  and  MAX  B  (fig.  8).  Operation  B  is 
as  ng  for  or  accepting  more  parts  than  A  is  producing  at  that  moment.  The  ver¬ 
tical  distance  on  the  chart  indicates  B  will  ask  for  21  more  parts  than  A  is 
producing  during  the  time  interval.  If  these  parts  are  to  be  provided  by  the  MHS 
so  as  not  to  adversely  affect  throughput,  then  21  parts  must  be  on  the  material 
handling  equipment  between  operations  A  and  B  to  buffer  these  operations.  There¬ 
fore,  the  material  handling  space  between  operations  A  and  B  should  be  20  vacan- 
cies  plus  21  fills  or  41  spaces. 


This  buffer  is  small  because  we  selected  two  operations  that  are  closelv 
balanced.  The  buffer  is  affected  according  to  the  unbalance  between  opera¬ 
tions.  Nevertheless,  the  methodology  indicates  the  buffer  required.  These  cal¬ 
culations  are  programmed  on  a  computer  and  the  charts  are  printed  out.  A  review 
of  these  charts  quickly  points  out  unbalanced  conditions  between  operations  and 
any  problems  that  may  result  in  buffer  sizes.  The  buffer  sizes  obtained  from 
this  methodology  can  then  be  used  as  input  for  buffer  sizes  to  the  GENMOD  program 
and  the  facility  tested  for  monthly  production  rates. 

Was  found  that  by  usinS  actual  calculated  buffer  sizes  as  input  to  the 
GENMOD  program,  more  realistic  results  were  immediately  obtained  based  on  what 
would  be  expected  from  past  experience  of  similar  lines  with  similar  quantities 
of  machine  tools.  Of  most  interest  was  the  finding  of  the  relationship  between 
u  er  size  and  predicted  monthly  production  capabilities  from  the  GENMOD  pro¬ 
gram.  A  buffer  size  of  one-hour  production  for  MSAAP,  260  parts  per  hour,  ap¬ 
pears  to  be  more  than  adequate  to  meet  most  conditions.  It  also  appears  that 

even  one-half  hour  production  would  be  adequate  in  most  cases  to  meet  monthly 
production  rates  provided  a  reasonable  maintenance  program  is  implemented  to  keep 
the  equipment  operating  at  reasonable  availability  values  and  catastrophic  down¬ 
time  was  kept  to  one  hour  supplemented  by  banks  in  the  line  as  indicated  bv 
GENMOD.  Dy 


results  of  the  GENMOD  also  indicated  where  cost  trade-off  might  be  made 
the  MHS  buffer  and  an  additional  machine  tool  to  better  balance  the  line 
and  this  factor  will  also  be  considered  in  future  design  of  metal  parts  produc¬ 
tion  line. 


CONCLUSIONS 


i.  It  appears  that  metal  parts  production  facilities  for  ammunition  projectiles 
can  be  simulated  by  computer  programs.  The  GENMOD  program  developed  at  ARRADCOM 
appears  to  be  the  most  suitable  for  such  simulations. 
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2.  Effective  simulation  programs  require  fairly  accurate  input  data  with  regard 
to  mean  time  between  failure  (MTBF),  mean  time  to  repair  (MTTR) ,  and  buffer  size. 

3.  When  provided  with  good  input  data,  the  GENMOD  program  appears  to  predict 
monthly  production  rates  that  coincide  with  past  experience  for  similar  produc- 
tion  lines. 


4.  Simulation  programs  such  as  GENMOD  permit  line  and  sensitivity  analysis  of 
parameters  previously  not  available  to  the  productin  line  designer.  Such  analy¬ 
sis  permits  the  designer  to  optimize  the  production  line  for  throughput  and  fa- 
cility  cost  before  it  is  constructed. 

5.  Line  simulation  permits  analysis  of  cost  trade-offs  for  integrated  facilities 
w  ere  material  handling  systems,  production  equipment,  etc.,  can  be  traded  off 

one  for  the  other  for  optimum  cost  while  still  maintaining  monthly  production 
r  a  t  e  s . 


6.  Simulation  programs  such  as  GENMOD  point  out  the  importance  of  preventive 
maintenance  programs  by  predicting  the  effect  on  monthly  production  caused  by 
catastrophic  failures.  The  trade-off  of  additional  equipment  and  storage  of 
parts  through  banks  becomes  evident  with  the  analysis  of  a  simulation  printout. 

7.  Analysis  of  GENMOD  printouts  permits  the  designer  to  assess  choke  points  in 
the  production  line  through  the  filling  or  emptying  of  buffers.  The  time  to  fill 
or  empty  such  buffers  can  also  be  seen  and  evaluated. 

8.  The  methodology  developed  using  the  binomial  distribution  is  very  helpful  in 
calculating  accurate  buffer  requirements. 

9.  The  binomial  distribution  methodology  also  illustrates  the  advantage  of  sev¬ 
eral  machines  producing  at  a  low  rate  rather  than  a  few  machines  producing  at  a 

high  rate  with  regard  to  reliability  and  meeting  expected  monthly  production 
r  at  e  s . 


10.  The  binomial  distribution  methodology  is  based  upon  availability  data  which 
is  readily  obtained  in  the  metal  parts  equipment  industry.  There  is  agreement 
within  the  industry  as  to  the  reliability  of  such  data. 


RECOMMENDATIONS 


1.  MTTR  and  MTBF  continue  to  be  collected  on  all  metal  parts  production  equip¬ 
ment  to  improve  the  reliability  of  predictions  developed  under  the  simulation 
programs. 

2.  Simulation  programs  such  as  GENMOD  continue  to  be  studied,  evaluated,  and 
used  on  existing  and  future  metal  parts  production  lines  to  improve  the  perform¬ 
ance  of  the  simulation  program  and  to  gain  a  deeper  insight  into  the  performane 
of  the  production  line. 
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3.  The  GENMOD  program  should  be  used  in  all  future  production 
fore  the  lines  are  built. 

4.  The  binomial  distribution  methodology  developed  under  this 
available  to  the  material  handing  industry  for  future  use. 


line  layouts  be- 
project  be  made 
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Table  2.  105-mm  downtime  record  sheet 
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Table  5.  Daily  statistics  for  finish  turn  band-1  (Station  105  at  Norris) 
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NOTE:  465  parts/5707  min  =  4.8  parts  hr 


Table  6.  Daily  statistics  for  finish  turn  band-2  (Station  106  at  Norris) 
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NOTE:  4286  parts/7756  min  =  33.15  parts/hr 

Operation  Finish  Turn  Band  (305)  *  33.15  +  4.8  =  37.95  parts/hr 


Table  7.  Station  summary  for  Norris  Industries 
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WELD  BAND  -  3  125.6  23.1  0.79940 
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MTBF  MTTR  bility 

Operation  by  machine  (min)  (min)  (%) 
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Table  7.  (cont) 
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Table  8.  System  summary 
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Table  9.  Norris  industrial  unlimited  buffer  summary 
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Table  10.  Control  point  data 
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By  rounding  off  partially  made  parts  per  BTI,  numbers  may  be  slightly  over  maximum  limits  stated. 


Table  12.  Mississippi  Army  Ammunition  Plant  line  simulation  data 
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Table  12.  (cont) 
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Table  12.  (cont) 
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Table  12,  (cont) 
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Table  14.  Simulation  Summary  for  Mississippi  Army  Ammunition  Plant: 

Total  Total  Total  Average  Average 

No,  of  failures  downtime  uptime  availability  stop  time 

Category  machines  (min)  (min)  (min)  (%)  (rain) 
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Table  14.  (cont) 
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Table  15.  Mississippi  Army  Ammunition  Plant  buffer  analysis 
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Table  16.  Louisiana  Army  Ammunition  Plant  metal  parts  system  summary 
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Table  16.  (cont) 
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Table  17.  Louisiana  Army  Ammunition  Plant  unlimited  buffer  summary 
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Table  21.  Percentage  of  time  various  numbers  of  machines  are  operating  (operation  B) 
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M  =  np  =  ( 6 ) ( 0 . 8 )  =  4.8  machines 

(4 . 8 ) ( 27 . 0833)  =  130  parts  per  hour  throughput 

Design  production  rate  =  27.0833  parts  per  hour  per  machine 
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Table  22.  Binomial  distribution-operation 
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Table  23.  Binomial  distribution-operation 
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TIME  INTERVALS  (MINUTES) 
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BUFFER  SPACE  BETWEEN  EACH  OPERATION 
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Figure  6.  Binomial  distribution 
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C/ML/L  A  77  VE  EL  AAS<££>  77 ME  —  M/A/U  T£S 

Figure  7.  Cumulative  parts  produced  versus  cumulative  elapsed  time  (MAX  A  -  MIN  B) 
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CO  MV  LA  77  VE  EL  APSE&  T/ME  —  M/A/O  TEE 

Figure  8.  Cumulative  parts  produced  versus  cumulative  elapsed  time  (MIN  A  -  MAX  B) 
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